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We consider extensions of the Dirac equation with mass terms mi + i7 5 m.2 and imi + 7 5 ni2- The 
corresponding Hamiltonians are Hermitian and pseudo-Hermitian ("7 s Hermitian"), respectively. 
The fundamental spinor solutions for all generalized Dirac equations are found in the helicity basis 
and brought into concise analytic form. We postulate that the time-ordered product of field operators 
should yield the Feynman propagator (ie prescription), and we also postulate that the tardyonic as 
well as tachyonic Dirac equations should have a smooth massless limit. These postulates lead to sum 
rules that connect the form of the fundamental field anticommutators with the tensor sums of the 
fundamental plane-wave eigenspinors and the projectors over positive-energy and negative-energy 
£N| , states. In the massless case, the sum rules are fulfilled by two egregiously simple, distinguished 

functional forms. The first sum rule remains valid in the case of a tardyonic theory and leads to 
' the canonical massive Dirac field. The second sum rule is valid for a tachyonic mass term and leads 

O^l \ to a natural suppression of the right-handed helicity states for tachyonic particles, and left-handed 

q ■ helicity states for tachyonic spin-1/2 antiparticles. When applied to neutrinos, the theory contains a 

' free tachyonic mass parameter. Tachyons are known to be repulsed by gravity. We discuss a possible 

f \ ' role of a tachyonic neutrino as a contribution to the accelerated expansion of the Universe ( "dark 

energy" ) . 
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I. INTRODUCTION 



A. Generalized Dirac Equations: Mass Terms and Dispersion Relations 

Dirac is often quoted as saying in some of his talks that the equation that carries his name [1, 2] is "more intelligent 
than its inventor". Of course, it needs to be added that it was Dirac himself who found most of the additional insight. 
Here, we are concerned with of extensions of the Dirac equation which contain both tardyonic and tachyonic mass 
terms. Tardyonic (subluminal) mass terms lead to dispersion relations of the form E = \/p 2 + m 2 , whereas tachyonic 
mass terms lead to superluminal dispersion relations of the form E = \/p 2 — m 2 , where E is the energy and p is the 
momentum. The generalized, matrix- valued mass term M enters the Dirac equation in the form (ij^d^ — M) tp(x) = 0. 
The 7 M arc 4x4 matrices that fulfill the relations {7^, 7"} = 2<? MI/ where we choose the space-time metric as 
g Ml/ = diag(l, — 1, — 1, — 1). The denote the partial derivative d/dx^ with respect to the space-time coordinate 
x^ = (t,x). It is quite surprising that a systematic presentation of the solutions of the generalized Dirac equations 
(i7 M 9 M — M) ip{x) = 0, in the helicity basis [3], has not been recorded in the literature to the best of our knowledge. 
While the following discussion is somewhat technical, we believe that it will be beneficial to give their explicit form, 
in order to fix ideas for the following discussion. 

We set H = c = eo = 1 and use the Dirac matrices in the standard representation 
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and define a = 7°7*. For the ordinary Dirac theory, one has M = mi (one should say more precisely M = mi 3.4x4) 
with a real mass mi, 

(i7"0„-mi) V(z) = 0. (1.2) 



The dispersion relation is E = y/p 2 + m\. The corresponding Dirac Hamiltonian reads 

H {1) = a -p + Pmi . (1.3) 

Extensions of the Dirac equation with pseudoscalar mass terms that contain the fifth current have been introduced in 
the literature. In Ref. [4], it is shown that for a mass term of the form M = mi + i7 5 m2, the fermion propagator may 
obtain nontrivial gradient corrections already at the first order in derivative expansion, for a position-dependent mass. 
In that case, the fermion self-energy may contribute to a conceivable explanation for CP-violation during clcctroweak 
baryogenesis, as pointed out in Ref. [4] . We thus study the following, generalized form of the tardyonic (subluminal) 
Dirac equation, 

(i7" - mi - i 7 5 ma) %j){x) = . (1.4) 
The dispersion relation is E = \J p 2 + m\ + m\. The Hermitian tardyonic Hamiltonian operator reads as 

#(*) = a -p + Pmx + i f3 m 2 ■ (1.5) 

We may indicate a further motivation for our study; namely, the unitarity of the S matrix implies the existence of 
useful relations [5] for the even powers (m 2 ) 2n obtained upon expanding a one- loop amplitude, formulated with a 
mass term mi + i 7 s m 2 , in powers of m 2 . This implies that a better understanding of the tardyonic equation with 
two mass terms could be of much more general interest. 

It has not escaped our attention that the chiral transformation 

fiexp(ij 5 9) = ^cos 9 + //sin 9 — mi + i7 5 m 2 (1.6) 

connects the two Hamiltonians and H^' for mi = [i cos 9 and m 2 = [i sin 9, but it is computationally easier 
and more instructive to consider the real and imaginary parts of the mass term separately. 

Within a systematic approach to generalized Dirac equations with pseudoscalar mass terms, we also consider 
tachyonic (superluminal) mass terms of the form M = 7 s m which induce a superluminal dispersion relation E = 
y/p 2 — m 2 . The corresponding generalized Dirac equation has been named the "tachyonic Dirac equation" and reads 
as follow [6-10], 

(i7 M 0„ ~ 7 5 m) i/}(x) = . (1.7) 
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The corresponding Hamiltonian reads 



H5 = a ■ p + (3 7' 5 m . 



(1.8) 



The relation H§ = 7°ii 5 l ~(— r ) 7 has been given in Re Is. [9, 10]. However, it is much more instructive to observe 
that H§ is 7 5 Hermitian, i.e., H$ = j 5 H^j 5 . The concept of 7 s Hermiticity is known in lattice theory [11, 12] and is 
otherwise called pseudo- Hermiticity [13-23]. 

An obvious generalization of the tachyonic case contains an imaginary mass and a 7 5 mass term, 



and is 7 s Hermitian, H' = j 5 H' + j 5 . For 7712 = 0, Eq. (1.9) has been discussed in Refs. [24, 25]. 

It is our goal here to present the fundamental eigenspinors corresponding to the plane-wave solution of the equa- 
tions (1-2), (1.4), (1.7), and (1.9) in a unified and systematic manner. Furthermore, we discuss the second-quantized 
versions of the fcrmionic theories described by the generalized Dirac equations. Anticipating part of the results, we 
may point out that the massless Dirac equation "interpolates" between the tardyonic equations (1.2) and (1.4) and 
the tachyonic equations (1.7) and (1.9). For zero mass, helicity and chirality arc equal. Helicity and chirality "depart" 
from each other in very specific directions, when the tardyonic and tachyonic mass terms are "switched on" , as we 
shall discuss in the following. 



The tachyonic generalized Dirac equations (1.7) and (1.9) describe the motion of superluminal particles, which may 
either be important for astrophysical studies (neutrinos) or for artificially generated environments such as honeycomb 
photonic lattices in which pertinent dispersion relations become practically important [26]. Contrary to other, some- 
what catchy, claims, the existence of superluminal particles would not falsify Einstein's theory of special relativity [27], 
which according to common wisdom is based on the following postulates: (i) The principle of relativity states that 
the laws of physics are the same for all observers in uniform motion relative to one another, (ii) The speed of light in 
a vacuum is the same for all observers, regardless of their relative motion or of the motion of the source of the light. 
Predictions of relativity theory regarding the relativity of simultaneity, time dilation and length contraction would 
not change if superluminal particles did exist. Furthermore, as shown by Sudarshan et at (Rcfs. [28-31]) and Feinberg 
(Refs. [32, 33]), the existence of tachyons, which are superluminal particles fulfilling a Lorentz-invariant dispersion 
relation E 2 = p 2 — m 2 , is fully compatible with special relativity and Lorentz invariance. According to special rela- 
tivity, it is forbidden to accelerate a particle "through" the light barrier (because E = m/y/l — v 2 —> 00 for v — > 1), 
but a genuinely superluminal particle remains superluminal upon Lorentz transformation. Significant problems are 
encountered when one attempts to quantize the tachyonic theories, but again, as shown in Ref. [9], these problems 
may not be as serious as previously thought. In particular, the so-called rcintcrpretation of solutions propagating into 
the past according to the Feynman prescription [31] is a cornerstone of modern field theory. Furthermore, it has been 
shown in Ref. [9] that tachyonic particles can be localized, and equal-time anticommutators of the spin-1/2 tachyonic 
field involve an unfiltered Dirac-<5 [see Eq. (37) of Ref. [9]]. 

Despite these arguments, we can say that, from the point of view of fundamental symmetries, accepting a su- 
perluminal neutrino would be equivalent to an "an ugly duckling" . Adding to the difficulties, we notice that recent 
experimental claims regarding the conceivable observation of highly superluminal neutrinos have turned out to be 
false. One may point out that a relative deviation v = c (1+6) with S ~ 10~ 5 at E ~ 30 GcV, as claimed by some recent 
experimental collaborations, would correspond to a negative neutrino mass square in the order of ~ — (100 McV) 2 , if 
one assumes a Lorentz-invariant dispersion relation [34]. Still, there is at present no conclusive answer regarding the 
conceivable supcrluminality of at least one neutrino flavor [35-39], and it is intriguing that all available direct mea- 
surements of the neutrino mass square have resulted in negative expectation values, still compatible with zero within 
experimental uncertainty, whereas published experimental best estimates for the neutrino speed [40-43] have been 
superluminal, again still compatible with the speed of light within experimental error. The recent ICARUS result [43] 
is consistent with this trend [44-50]; the best estimate for the neutrino velocity is superluminal, but the deviation 
from v v = c is statistically insignificant. The OPERA collaboration [51] has indicated a preliminary, revised result of 
(v — c)/c~ (2.7 ± 3.1(stat)^2 g( s y s -)) x 10~ 6 - Neither subluminal nor superluminal propagation velocities are excluded 
based on the available experimental data. The "ugly duckling of a superluminal neutrino" is not beautiful; if we are 




(1.9) 



The dispersion relation is E = \Jp 2 — m\ — m\. The corresponding Hamiltonian reads 

H' = a ■ p + i /3 mi + /3 7 s 7712 , 



(1.10) 



B. Tachyonic Dirac Equation and Neutrinos: Possible Connections 



4 



to consider accepting it, then we should be able to hope that the emergence of at least one "swan" (or "intellectual 
benefit" ) should be the result of this operation. 

Before we discuss the possible emergence of these benefits, let us include some historic remarks. According to 
reliable sources [52], Professor J. A. Wheeler, in his later years at the University of Austin (Texas), used to argue 
that the neutrino has to be massless, necessarily, and that in his opinion, it could only be a massless Weyl particle 
with definite helicity (and chirality). Recently presented arguments [53] regarding the possibility of overtaking a 
subluminal, left-handed neutrino, looking back and seeing a right-handed neutrino, were supposedly already used 
by Wheeler in order to dispel the conceivable existence of a neutrino mass term. This paradox has been termed 
"autobahn paradox" in Ref. [53] and excludes a Dirac neutrino unless one assumes exotic processes like sterile right- 
handed massive neutrinos. [The problem with a right-handed sterile massive neutrino is that for massive neutrinos, 
chirality and helicity are different, hence a V — A coupling of the form 7 A '(1 — 7 5 ) no longer vanishes for massive Dirac 
neutrinos if one uses the canonical eigenstates of the massive Dirac equation. One therefore has to invoke additional 
exotic mechanisms in order to ensure the "sterility" of the right-handed Dirac neutrinos.] Wheeler also disliked [52] the 
notion of a Majorana neutrino, arguing that the charge conjugation invariance condition imposed on the Majorana 
particle precludes the existence of plane-wave solutions to the Majorana equation, and maximally violates lepton 
number. Again, these arguments [54] are in full agreement with those recently given in Ref. [53]. 

The original standard model thus called for manifestly massless neutrinos. The commonly accepted observation 
of neutrino oscillations precludes the possibility that all three generations of neutrino mass eigenstates are massless. 
Lepton number conservation is based on the global gauge symmetry tp — > -0 exp(iA), applied simultaneously to 
all lepton fields. A Majorana neutrino would destroy lepton number as a global symmetry but solve the "autobahn 
paradox" , because a Majorana neutrino would be equal to its own antiparticlc and thus, looking back, the right-handed 
neutrino state would consist of the same particle=antiparticle. 

On the other hand, if we assume that the neutrino is described by the tachyonic Dirac equation, then the following 
statements arc valid: 

• Statement #1: We can properly assign lepton number and use plane- wave eigenstates for incoming and outgoing 
particles, while allowing for nonvanishing mass terms and thus, mass square differences among the neutrino mass 
(not flavor) eigenstates. 

• Statement #2: There is a natural resolution for the "autobahn paradox" because a left-handed spacclikc neutrino 
always remains spacelike upon Lorentz transformation and cannot be overtaken. 

• Statement #3: The right-handed particle and left-handed antiparticle states are suppressed due to negative 
Fock-space norm. 

• Statement #4: At least qualitatively, tachyonic neutrinos could yield an explanation for a repulsive force on 
intergalactic distance scales as they are repulsed, like all tachyons, by gravitational interactions ("dark energy"). 

Pauli [55] postulated the existence of neutrinos, on the basis of the conservation of angular momentum and energy, 
and also introduced pseudo-Hermitian operators [13]. Here, we describe conceivable connections of neutrino physics 
and pseudo-Hermitian operators. Final clarification can only come from experiment. When in 1956, F. Reines and 
C. Cowan [56] discovered the electron neutrino, two and a half years before Pauli's death, Pauli replied [57] by telegram: 
"Thanks for message. Everything comes to him who knows how to wait. Pauli." In defense of the tachyonic hypothesis, 
we would like to stress that a tachyonic Dirac neutrino would allow us to retain lepton number conservation as a 
symmetry of nature. We would thus like to write up these thoughts in the current article, with attention to detail. We 
should point out that our approach fully conserves Lorentz invariance, in contrast to the extensions of the Standard 
Model based on Lorentz-violating terms which can otherwise lead to superluminal propagation (see Tables 11 and 13 
of Ref. [58]). 

Units with h = c = e = 1 are used throughout the paper. The organization is as follows: In Sec. II, we discuss 
massless and tardyonic theories. Tachyonic extensions of the Dirac equation are discussed in Sec. III. Connections of 
the fundamental tensor sums over the eigenspinors with the derivation of the time-ordered propagator are analyzed 
in Sec. IV. A candidate for dark energy is presented in Sec. V. Conclusions are reserved for Sec. VI. 

II. GENERALIZED DIRAC EQUATIONS: MASSLESS AND TARDYONIC THEORIES 

A. Massless Dirac Theory 



The massless Dirac equation and the massless Dirac Hamiltonian read as 

vy^d^ ip (x) = 0, H = a-p. 



(2.1) 
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We note that Hq is both Hermitian as well as 7 5 Hermitian, i.e., Hq — 7 s Hq j 5 . The dispersion relation is E = \k\. 
With k^ = (E, k), we seek positive-energy and negative-energy solutions of the form 



tp(x) = u a {k) exp(— ik ■ x) , <p(x) — v a (k) exp(i/c • x) , 



(2.2) 



where a = ± denotes a quantum number which is equal to the helicity for positive-energy states, and equal to the 
negative of the helicity for negative-energy states. With ji = 7 M /c M , we have ^u±{k) = ftv±(k) = 0. In the massless 
limit, the solutions to the Dirac equation are given as (see Chap. 2 of Ref. [59]) 



u + (k) 



_l_ { a + (k) 
y/2 \ a+(k) 

v+(k) = — u + (k) , v-(k) 

The well-known helicity spinors are recalled as 

a+(k) -- 



, u_(fc) 




-U-(k) 



a-(k) 



sin (|) 



(2.3a) 
(2.3b) 

(2.4) 



These fulfill the fundamental relations (er ■ k) a a (k) — aa a (k), as well as ^2 a <y{k) (£) a+ (k) = 3.2x2 , an d J2a a a a{k)(S> 

a+(/c) = a ■ k, where k = k/\k\ and the sum over a is over the values ±1. The sums over the fundamental bispinors u 
and v fulfill the following sum rules, 



rule I: 



2\k\u a {k) ®u a {k) Y 2|fcK(fc)<g)u CT (fc) 



as well as 



Sum rule II: 



2\k\{-<y)u ry {k)®u ty {k) 1 b J2 2|fc|(- ( r)t; (T (fc)®TJ CT (fc)7 5 = ^. 



(2.5) 



(2.6) 



Sum rule I can be obtained by a quick explicit calculation, and sum rule II holds because in the massless limit, helicity 
equals ichirality (positive sign for positive energy, negative sign for negative energy). We denote the Dirac adjoint as 

u a (k) — w+(fc)7°. One can easily check by an explicit calculation that u a (k)^f 5 = (^y 5 j° u a (k)j = (—a)u a {k) and 

TJ CT (/c)7 5 = ^7 5 7 v a (k)J — (—<7)v a (k). We can thus introduce a factor (— a) 2 — 1 under the summation over spins in 
Eq. (2.5) and replace one of the factors (—a) by a multiplication of the Dirac adjoint spinor from the right by the fifth 
current. The Lorentz-invariant normalization of the massless solutions vanishes, i.e., u a (k) u a (k) — Vcr(k) v a (k) = 0. 

Eigenstates of the massless Hamiltonian Hq — a ■ p have to be eigenstates of the chirality operator 7 5 because the 
chirality commutes with the Hamiltonian, in the sense that [7 5 ,i?o] = 0. Furthermore, 



H = a ■ p = |p| 7 5 



(2.7) 



where S = 7 s a is the vector of 4 x 4 spin matrices, and the helicity operator is identified as E • p/\p\. Let Ai be the 
eigenvalue of chirality and A2 be the eigenvalue of the helicity operator. Then, the eigenvalue of the Hamiltonian is 
Eq = \p\ Ai A2. Since Ai = ±1 and A2 = ±1, we easily recover the known fact that helicity equals chirality for positive 
energy, whereas the relation is reversed for negative-energy states (see also Chap. 2.4 of Ref. [59]). We are aware of 
the fact that the considerations reported in the current section partly refer to the literature but we give them in some 
detail because they are essential for the following considerations. 

B. Massive Dirac Theory 

We start from the ordinary Dirac equation given in Eq. (1.2), which reads (i7 M <9 M — mi) iJj(x) — 0. In the helicity 
basis, the fundamental spinor solutions read as 



ip(x) = U±(k) exp(— ik ■ x) , <f)(x) = (k) exp(ifc • x) . 



(2.8) 
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The algebraic relations that have to be fulfilled by the bispinor amplitudes U± (k) and V±' (k) read as follows 



-mi) U£\%) =0, (ft + mi) V± L) (k) = 



(2.9) 



The dispersion relation is E^ = y k 2 + m\. In the helicity basis, the solutions of the equation (2.9) with a tardyonic 
m\ mass term are easily written down, using the identity (ft — mi)(ft + mi) = k 2 — m\ = 0. With an appropriate 



normalization factor, and after some algebraic simplification, the positive-energy solutions read as follows, 

(ft + mi) u+(k) 



U?\k) 



U W (k) = 



y/(EW - \k\y+m 2 
(^ + mi)ii,(l) 

(J5(l) _ |fc|)2 +m 2 



VV 2£?(D a+[k) ) 



V 



2E^ 



a_ (k) 



mi 



2EW 



a-(k) 



(2.10a) 



(2.10b) 



The negative-energy eigenstates of the tardyonic equations in the helicity basis are given as 



vi 1] (k) 



vP(k) = 



(mi - ft) v+(k) 
(EW - \k\) 2 + mf 

(mi — ft) v-(k) 
( E W - \k\) 2 +m\ 



a + (k) 



\ 



EV 



2EW 



a + {k) 



( M 



2EW 



1711 ( M \ 



V 



E^ 



mi 



2EW 



a_(fc) 



(2.11a) 



(2.11b) 



These solutions are consistent with those given in Ref. [3] and in Chap. 23 of Ref. [60], and the normalizations are 
(a = ±) 



uP+Qs) uP(k) = VP+(k) vP(k) = 1 . 
One can change the normalization according to 

u£\k) = ( ) up® , VP® = ( — ) V^(k) . 



\ mi 



V m i 



(2.12) 



(2.13) 



The Lorentz-invariant normalization is equal to one for the fundamental positive-energy bispinors and equal to minus 
one for the fundamental negative-energy bispinors, 



u¥(k)uP(k) = i, V ( a L> (k)vP(k) 



-1 



A little algebra is sufficient to reproduce the following known sums over bispinors, 



ft + mi 
2mi 



J2vP(k)®v ( ^(k) 



ft — mi 
2mi 



(2.14) 



(2.15a) 



In accordance with general wisdom about the tardyonic case, these do not involve a helicity-dependent prefactor. The 
sum rule (2.15) is of type I [see Eq. (2.5)]. 
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C. Two Tardyonic Mass Terms 



Inspired by the discussion in Sec. I, we consider an equation with two tardyonic mass terms, which has already 
been indicated in Eq. (1.2) and reads (iry^d^ — mi — 1^1712) tp^(x) — 0. For the corresponding bispinors in the 
fundamental plane-wave solutions, this implies that 



Qi - mi - i7 5 m 2 ) U^\k) = , 
-mi- i7 5 m 2 ) V^\k) = . 



(2.16a) 
(2.16b) 



The dispersion relation is E^> = \ k 2 



m\. The fundamental positive-energy bispinors read as follows, 



Uf{k) 



+ mi — i7 5 7712) w+(fc) 



(EW - \k\) 2 +mj+m 



( mi - im 2 + £ (t) - \k\ a+(k) \ 

(£(*)- \k\Y + m\ + m\ ^ 
mi — i 77i2 — E^' + \k\ a+(k) 

( mi + im 2 + E® - \k\ a_(fc) \ 



\ J(EW - |fc|) 2 +777,2 +m 2 



(2.17a) 



^(jfc) 



(ft + mi — i7 5 m 2 ) U- (fc) 



(J5W - |fc|) 2 +mj +m\ 



V2 



m - im 2 + E®> - |fc| a_(Jfe) 



V2 



'(£(*) - |jfc|) 2 +7nf + m 2 

\^{EW -\k\y+ml + l „ 2 
The negative-energy eigenstates of the equation with two tardyonic mass terms are given as 



(2.17b) 



/ -mi + im 2 + E^ -\k\ a+(k) \ 



T/ (t),A (-/^- i7 5 77i 2 +mi)u + (fc) 



(EW - \k\) 2 + m\ + m 



(EW - \k\) 2 +m 2 +m 2 



V2 



-mi + im 2 - E^> + \k\ a+(k) 



7 2 V2 



\ y/(EW - \k\) 2 +m 2 +m 
for negative helicity (positive chirality in the massless limit) and 

/ -mi-im 2 + E& -\k\ a_(fc) \ 



r (i)/A _ (-<^-i7 5 "72 +mi)^_(fc) 



Vi ; (fc) = 



(£« - |fc|) 2 +m 2 + m 2 



(SW- |fc|) 2 + mf + ml 



V2 



mi + im 2 + E^ - \k\ a_(fc) 
_ __ 



\ J{EW -\k\) 2 +m 2 + ' 



) 



for positive helicity (negative chirality in the massless limit). In the massless limit (first E® — > \k\, then 7772 — > 0, and 
then mi — > 0), we again reproduce the massless solutions, U+(k) — > u + {k), U_(k) — > u_(fc), V+Qc) — > v + (k) and 

V W {k) — > V-(k). Of course, in the limit mi — > 0, one also has to expand the normalization denominators in powers of 
mi. For 7772 = the solutions (2.17) and (2.18) reduce to the solutions of the ordinary Dirac equation in Eqs. (2.10) 
and (2.11), and can be expanded in 7772 to yield corrections to the ordinary Dirac equation for small m 2 , i.e mi 3> m 2 . 
The states are normalized with respect to the condition 



U^(k)U^(k) = V^(k)V^(k) = l 



In the normalization 



U$\%) = 



1/2 



£#>(*), V«(fc) = 



1/2 



V^(k), 



(2.19) 



(2.20a) 



the positive-energy solutions acquire a "positive Lorentz-invariant norm" , whereas the negative-energy solutions have 
"negative Lorentz-invariant norm" , 



-1 



(2.21) 
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After some algebra, one can derive the following sums over bispinors, 



2toi 



ft — mi + 17 5 m 2 
Inii 



(2.22) 



These are easily identified as the positive- and negative-energy projectors. The sum rules do not involve helicity- 
dcpcndcnt prcfactor and are of type: I [sec: Eq. (2.5)]. 

The solutions (2.17) and (2.18) approach the massless solutions if one replaces ni2 — > first and then lets nil — > 0. 
They are thus useful for systems where the nil mass is greater than 77x2- For 7712 3> mi, one would like to calculate 
solutions that approach the massless case for the sequence mi — > 0, then — > 0. These read as follows, 



U'W (k)=ia V® (k) , V™ (k) = ia V® (k) . 
In comparison to the solutions (2.17) and (2.18), they acquire a nontrivial phase factor. 



(2.23) 



III. GENERALIZED DIRAC EQUATIONS: TACHYONIC MASS TERMS 

A. Tachyonic Dirac Equation 

The tachyonic Dirac equation is given in Eq. (1.7) and reads (ij^ — 7 5 mj ip(x) = 0. The fundamental bispinors 
entering the equations fulfill the equations 



(ft-^m) U±(k) = 0, (ft + ~f 5 m) V±(k) =0. 



(3.1) 



Using (ft — 7 5 m) (ft — 7 15 m) = k 2 + m 2 and some algebra, the prefactors in the fundamental bispinors (for positive 
energy) take a very simple form, 



( m 



U+(k) = 



(7 s m — ft) u+(k) 
(E- \k\) 2 +m 2 

U (k)= (^7 5 m)._(fc) 
y/(E- \k\) 2 +m 2 

For negative energy, the solutions read as follows, 

(7 5 m + ft) v+(k) 



2\k\ 



m fu\ ^ 
a + (k) 



\k\- 



a + (k) 

V V 2|fc| v 



( ME™ a.® ^ 
V 2|fc| 

v V 2 * ' j 



V+(k) 



V-(k) = 



(E-\k\) 2 +m 2 



(-ft — 7 5 m) v-(k) 



(E-\k\) 2 +m 2 



l \k\-m - \ 

l^\kT a+{k) 

j \k\+m -> 



2|fc| 
Ifcl — m 



a-(k) 



V V 2|fe| ~ v J 



The normalization condition is (a = ±) 



(3.2a) 



(3.2b) 



(3.3a) 



(3.3b) 



U+(k)U a (k)^V+(k)V a (k) = 1. 



(3.4) 
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One can change to Lorentz- invariant normalization by a multiplication with (jfcl/m) 1 / 2 , 



U a (k) 



1*1 
m 



1/2 



1/2 



t/ CT (fc), V CT (fc) 



(3.5) 



The "calligraphic" spinors fulfill the following helicity-dependent normalizations, 

U a (k) U a {k) = a , VS) V«r(£) = -a , 



(3.6) 



where we observe that a is a good quantum number because the helicity operator commutes with the Hamiltonian (1.8). 
The sum rule fulfilled by the fundamental plane- wave spinors is of type II [see Eq. (2.6)] . For the positive-energy spinors, 
we have 





^(~<j)U a (k)®U a (k) 7 5 

a 


H 


— 7 5 TO 
2m ' 


where for the negative-energy spinors, the sum rule reads 






H 


■f 7 5 TO 

2m 



(3.7a) 



(3.7b) 



The expressions on the right-hand sides are the positive- and negative-energy projectors. 



B. Two Tachyonic Mass Terms 



We study the equation (i^f^ — im^ — 7 5 m 2 ) ip(x) = 0, as given in Eq. (1.9). The fundamental spinors, which we 
denote as U'±(k) and V±(fc), fulfill the following equations, 



(^-im 1 - 7 5 m 2 ) U^(k)=0, 
(/f + im 1 + 7 5 m 2 ) Vl(k) = 0. 

The positive-energy solutions are obtained using the identity (ft — imi — 7 5 m 2 ) (ft + imi — 7 5 m. 2 ) 
With E' = \/k 2 — m? — mf, they read as follows, 



U' + (k) = 



U'_(k) = 



/ 


imi 4 


- m 2 - 


£'4 


\k\ 


a+(k) \ 




vV- 


\k\r 


4- m\ 


+ m\ 






imi 4 


■ m 2 -+ 


-E'- 


\k\ 


a+(k) 


V 




\k\) 2 


4- m\ 


+ m\ 


V2 J 




imi ■+ 


m 2 4 


E' - 


\k\ 


a_(k) \ 






\k\) 2 


4- m\ 


+ m| 






—imi - 


- m 2 - 


IE' - 


-\k\ 


a-(k) 


V 


J{E'- 


\k\f 


+■ m\ 


+ m\ 


V2 J 



(3.8a) 
(3.8b) 



(3.9a) 



(3.9b) 
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The negative-energy solutions for the tachyonic equation with two mass terms, are given as 



vUk) 



( 


imi — 


m 2 — 


E' + 


\k\ 


0+(fc) \ 




M- 




f m\ 


+ m% 


V2 




imi — 


m 2 - 


E' - 


\k\ 


a+(k) 


V 1 


M- 


\k\f- 


f m\ 


+ m\ 


V2 j 


/ 


—imi " 


- m 2 - 


YE' - 


-\k\ 


a_(k) \ 




M- 


\k\y- 


f m\ 


+ ml 






imi 4 


m 2 + 


E' - 


|fc| 


a-(k) 


V 


M- 


\k\y- 


f m\ 




V2 J 



(3.10a) 



(3.10b) 



The normalization condition is U'+{k) U' a {k) = V ( ! r + (k) V^(fc) = 1. We use a definition of the "calligraphic" spinors 
analogous to Eq. (3.5), 



/ ,-, \ 1/2 / \ 1/2 

\ 7772 / \ 7772 / 



(3.11) 



In analogy to Eq. (3.7), a sum rule of type II [see Eq. (2.6)] is fulfilled by the fundamental plane-wave spinors, 



]r(-a)^(fc)^(fc) 7 5 = 



+ imi — 7 7Ti2 
2m 2 



5 ft — i 77X1 + 7 5 777 2 



2m 2 



(3.12) 



We thus obtain the desired projectors onto positive- and negative-energy solutions for the Dirac equation with two 
tachyonic mass terms (1.9). The generalized equation (ij^d^ — imi — l 5 m 2 s ) ip(x) — is fully compatible with the 
Clifford-algebra based approach recently described in Ref. [61]. 



IV. THEOREMS FOR GENERALIZED DIRAC FIELDS 



A. Spinor Sums and Time-Ordered Propagator 



Our central postulate regarding the quantized fermionic theory is that the time-ordered vacuum expectation value 
of the field operators should yield the time-ordered (Feynman) propagator, which, in the momentum representation, 
is equal to the inverse of the Hamiltonian (upon multiplication with 7 ). This postulate implies that under rather 
general assumptions regarding the mathematical form of the elementary field anticommutators, sum rules have to be 
fulfilled by the tensor sums over the fundamental spinor solutions. It is perhaps not surprising that these sum rules, 
are precisely of the form investigated in Sees. II and III of this paper. 

For definiteness, we consider the solution of the tachyonic Dirac equation (Sec. Ill A). The generalization to other 
generalized Dirac equations is straightforward. We start from the field operator [9] 

/ 1 E {K{k)U a {k)e-^ + dt{k)V,{k)e^} , 

k = (E,k), E = Ej: = \Jk 2 - m 2 - i e , (4.1) 
where b a annihilates particles and <f+ creates antiparticles. The following anticommutators vanish, 



{b a (k),b p (k')} = {b+(k),b+(k')}=0, 
{d a (k),d p (k')}= {d+{k),d+(k')} =0, 



(4.2a) 
(4.2b) 



11 



We assume the following general form for the nonvanishing anticommutators, 

{b a {k),b+{k')} =f(a,k) (2tt) 3 E-S 3 (k-k')S ap , (4.3a) 

{d a (k) ) d+(k')}=g(a,k) (2iT) 3 -S 3 (k~k')5 ap , (4.3b) 

m 

with arbitrary / and g functions of the quantum numbers a and k. One might argue that since / and g must be 
dimensionless, they can depend only on the dimensionless arguments a and k/m, but that is a detail of the discussion 
which we do not pursue any further. Our only assumption concerns the fact that the field anticommutators should be 
diagonal in the helicity and wave vector quantum numbers, leading to the corresponding Kronecker and Dirac-<5's. 

We assume that the spin-matrix T either constitutes a Lorentz scalar or a pseudo-scalar quantity, which is a scalar 
under the proper orthochronous Lorentz group. The time-ordered product of field operators reads as 

(Q|T#r)^(y)r|0)= J {&(x° - y°) e - ifc '^) £ f(a, k) U a (k) ® U a {k) T (4.4) 



0(y° - x°) e ik <*-^ Yl 9(<r, k) V a (k) ® V a (k)T} 



This equation contains the same coefficient functions / and g that enter into Eq. (4.3). In order to proceed with the 
derivation of the propagator, we must postulate that the following sum rules hold, 

/fo fc) U a {k) ®U a {k)T= ^ ~^ m , ]T g(a, k) V a (k) ® V g (fc) V = * + ^ m . (4.5) 

a a 

The sum rule (4.5) is crucial for the further steps in the derivation of the time-ordered propagator. Introducing a 
suitable complex integral representation for the step function, one obtains from Eq. (4.4), using Eq. (4.5), after a few 
steps which we do not discuss in further detail [see also Eqs. (3.169) and Eq. (3.170) of Ref. [59]] 

(0 \T^(xU(y)T\ 0) =i / — ™ [ ^ e -^( X "-y")+iHS-y) 1^-^-k-^m 
\ U I WWJ-H u ; l ] (27r) 3£ J 2ir e 2m(k -E + ie) 

d3fc m f dk -ik -(x -y°)+ik-(x-y) ~l° k + 7 \ k + "f 5 TU 

(2tt) 3 E J 2tt 2m(fc + E - ie) ' K ' ' 



The convention is that in any integrals J d 3 fc, the component fco is set equal to E — \/k 2 — m? in the integrand when 
it occurs in scalar products of the form k ■ (x — y) etc., but if the integral is over the full d k, then the integration 
interval is the full fco € (—00, 00). With the convention which is adopted in many quantum field theoretical textbooks, 
including Refs. [59, 62], we finally obtain the result 

x 1 J (27r) 4 k z + m £ + ie 



7 5 m 



The tachyonic propagator St is identified, under the integral sign, as 

Mfc) - s rr ; — t = ■ ( 4 - 8 ) 

w ^-7 5 (TO + ie) k 2 +m 2 + ie v ' 

The sum rule (3.7) implies that the derivation is valid for the choice 

f(a,k)=g(a,k) = -a, T = 7 5 , (4.9) 

in which case the relations given in Eq. (4.5) are fulfilled. Note that this observation does not imply that the choice (4.9) 
necessarily is the only one for which we are able to fulfill the postulates given in Eq. (4.5), but it is the only structurally 
simple choice that we have found. 

For the egregiously simple choice (4.9), let us study the transition to the massless limit (2.6) in some further detail. 
Indeed, in the limit m —¥ 0, the denominator of the spin sums in Eq. (3.7) vanishes, and a finite limit is obtained after 
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multiplication with 2m, 

lim V2m(-cr) W CT (fc) <g> Z7 ff (fc) 7 5 = fz, (4.10a) 
lim V 2m (-cr) V CT (fc) ® VJk) 7 5 = (4.10b) 

cr 

In order to compare the normalizations of the fundamental spinors in the massless limit, we calculate the following 
quantities, 

Z7 CT (fc) 7 5 7°W cr (fc) = -a-, u r7 {k) 1 5 1 °u rT (k) = -a, (4.11a) 

m 

V a {k) 1 b 1 °V a {k)= -a — , zJ (T (fc) 7 5 7 % cr (fc) = -a, (4.11b) 

m 

where the fundamental spinors u a and v a of the massless equation have been given in Sec. II A. Observing that 

E = \k\ in the massless limit, the identifications ^JmU a (k) — > \J\^\u a (k) and y/mV a (k) — > \J\^\v a (k) : as implied by 
Eq. (4.11), show that the identities (4.10) precisely reduce to the sum rule (2.6) in the massless limit. 

B. Generalized Field Anticommutators for Tardyonic and Tachyonic Fields 

For dcfinitcness, we have considered the case of the tachyonic Dirac field in the above derivation. The decisive 
observation is that the choice 

tachyonic choice: /(a, k) = g(a, k) = —a , T = 7 5 , (4-12) 

is consistent with both massive tachyonic fields discussed in Sees. Ill A and IIIB, whereas the tardyonic choice 

tardyonic choice: /(cr, k) = g(a, k) = 1 , T = 1.4x4 , (4-13) 

yields the time-ordered propagator for both massive tardyonic fields discussed in Sec. II B and II C. The nonvanishing 
anticommutators for tardyons take the simple form [cf. Eq. (4.3)], 

tardyonic anticommutators: 

{&„(*), &+(«/)} = (2nf-S 3 (k-k')5 ap , {d a (k),d+(k')} = (2t:) 3 -S 3 (k-k')S ap . (4.14) 
m m 

Again, compared with Eq. (4.3), the nonvanishing anticommutators for tachyons take the simple form 
tachyonic anticommutators: 

{Mfc),6+(fc')} = (-0 (2?r) 3 -S 3 (k - k')8 ap , |d a (fc),d+(fc')} = {-a) (2tt) 3 -§\k - k')S up . (4.15) 

m H m 

With these universal choices, the theory of the tardyonic and tachyonic spin-1/2 fields can be unified. The time-ordered 
propagator is given as 

(0 \T^x)ii{y)T\ 0) =iS(x-y), S(x ~v) = j ^ S(k) . (4.16) 

We use the sum rules for the tensor sums over fundamental spinors given in Eq. (2.15) (for the tardyonic Dirac 
field), Eq. (2.22) (for the tardyonic Dirac field with two mass terms), Eq. (3.7) (for the tachyonic Dirac field) and 
in Eq. (3.12) (for the Dirac field of imaginary mass). Going through the exact same derivation as outlined above 
in between Eqs. (4.4) and (4.8), we obtain the following results for the time-ordered propagators of tardyonic and 
tachyonic fields. For the tardyonic Dirac field (Sec. II B), one has 

W = g 1 ■ = J + T ■ (4-17) 
ft - mi + ie k 2 - m\ + l e 
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For the tardyonic field with two mass terms (Sec. (II C)), the Feynman propagator is easily found as 

= Tl = j + m r l7 2 5W2 , (4-18) 

ft — mi +ie — i 7° (m 2 — ir]) k z — m\ — m 2 + i e 

where e and r\ are infinitesimal imaginary parts. Both tardyonic mass terms acquire an infinitesimal negative imaginary 
part, and the prefactor m/E from Eq. (4.1) in the field operator needs to be replaced by m\/ E^> where the tardyonic 

energy is = \Jk 2 + m\ + m\ ■ For the tachyonic Dirac field (Sec. Ill A and Ref. [9]), one has the result given in 
Eq. (4.8). Finally, for the Dirac field with two tachyonic mass terms, we have 

S'(k) = / 2 +i 7~t m2 (4-19) 
w k 2 + m\ + m\ + 1 e y ' 

In the latter case, the prefactor m/E from Eq. (4.1) in the field operator needs to be replaced by mi/E' where the 



tachyonic energy is E 1 = y k 2 — m\ — m\ . For both tachyonic fields discussed here, the mass acquires an infinitesimal 
positive imaginary part, as manifest in the results given in Eqs. (4.8) and (4.19). 



C. Tachyonic Gordon Identities 

It is useful to illustrate the derivation outlined above by exploring its connection to tachyonic Gordon identities. 
For dcfinitcness, we again concentrate on the tachyonic Dirac equation discussed in Sec. Ill A. The matrix element of 
the vector current finds the following Gordon decomposition for positive-energy spinors, 

U ± {k')rU±(k') = -Lw±(fc')7 5 [(fc^-n+i^ {k' v + k v )]U±{k) . 
Zm 

For negative-energy solutions, the identity reads as 

1 



V ± (fc') 7 M V ± (fc) = --^V ± (fc')7 5 [(k'^-k^ + ia^ (K + K)] V±(fc). 
2m 



For k' = k, one has 



W±(Jb)yW±(jfe) = —U ± {k)^a^k v U ± {k), (4.21a) 



m 



V±(k)^V±{k) = - —V ± (k)^ a^k y V ± {k). (4.21b) 
m 



The matrix element of the axial current reads 



U±(k')j 5 ^U±(k) = -^-U±{k') [(k'» + k») + ia» v {k' v - K)]U ± {k) , (4.22a) 
2m 



whereas for negative-energy solutions 



V±(fc') 7 VV±(fc) = ^-V±{k')[(k' fi + k>*)+i ( T^{k' u -K)}V±{k). (4.22b) 



For k! = k, this simplifies to 

U± (fc) 7 5 7 ' 1 W± (k)= - —U± (k) k" U±(k), (4.23a) 
m 

V± (fc) 7 5 7 AI V± (k) = — V± (k) k" V± (k) . (4.23b) 
m 

The results (4.22) and (4.23) for the tachyonic axial vector current have a similar structure as the Gordon decomposi- 
tion for the tardyonic vector current obtained with the ordinary Dirac equation [see Eq. (2.54) of Ref. [59]]. The role 
of the Dirac adjoint for the tardyonic case is taken over by the "chiral adjoint" U±(k)'y 5 for the tachyonic particle. 
Here, the designation "chiral adjoint" is inspired by the fact that U±(k)^f 5 U±(k) transforms as a pseudo-scalar under 
Lorentz transformations. 
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The structure of Eqs. (4.21) and (4.23) is somewhat peculiar with regard to parity. In Eq. (4.21a), an apparent 
vector current on the left-hand side appears to transform into an axial current on the right-hand side, whereas in 
Eq. (4.23a), an apparent axial vector on the left-hand side of the equation becomes what appears to be a vector on 
the right-hand side. The reason lies in the more complicated behavior of the tachyonic Dirac equation under parity 
as investigated in Ref. [10]. Namely, the tachyonic Dirac equation (1.7) contains a term which transforms as a scalar 
under parity, 

i l v ^ A 7° (i 7° d + i f (-9,)) 7° = i l" d„ (4.24a) 
as well as a term which transforms as a pseudoscalar, 

7 5 to -A 7° (7 5 to) 7° = - 7 5 m . (4.24b) 

The mass term in the tachyonic Dirac equation is pseudoscalar and changes sign under parity. Indeed, in Ref. [10], 
the tachyonic Dirac equation has been shown to be separately CV invariant, and T invariant, but not V invariant, 
due to the change in the mass term. 

In order to put this observation into perspective, we recall that the entries of the electromagnetic field strength 
tensor are composed of axial vector components (magnetic B field), as well as vector components (electric E field). 
The transformation properties of the electromagnetic field strength tensor under the proper orthochronous Lorentz 
group are nevertheless well-defined. 

The transformation (4.24b) can be interpreted as a transformation to — > —to under parity. Thus, if we interpret 
the mass m as a pseudoscalar quantity, then the right-hand sides of (4.21a) and (4.23a) transform as a vector and 
an axial vector, respectively. It is the parity non-invariance of the mass term in the tachyonic Dirac equation which 
leads to the somewhat peculiar structure of Eqs. (4.21a) and (4.23a). 



D. Helicity— Dependence and Gupta— Bleuler Condition 

The anticommutator relations for tardyons given in Eq. (4.14) imply that both left-handed as well as right-handed 
helicity states, for both particles as well as ant ip articles, have positive norm. However, the anticommutator relations 
for tachyons given in Eq. (4.15) imply that right-handed particle as well as left-handed antiparticle states acquire 
negative norm. This is shown in Eqs. (31) and (32) of Ref. [9]. Indeed, for one-particle states |lfc )CT ) = b+(k)\0), 

(U,Mk,a) = (0 \b a (k) b+(k)\0) = (0\{b a (k),b+(k)}\0) = (-a) V ~ , (4.25) 

where V = (27r) 3 <5 3 (0) is the normalization volume in coordinate space. The Fock-space norm (lk,a\^k,cr) is negative 
for cr = l. 

For clarification, the corresponding Gupta-Bleuler condition should be indicated explicitly. In full analogy to the 
instructive discussion of the Gupta-Bleuler mechanism for the photon field, as given in full clarity in Chap. 9b of 
Schweber's textbook [63], we select the positive- and negative-frequency component of the (er = l)-component of the 
neutrino field operator, 

= / | W U*=i&) e~ ik x , (4.26a) 

^= l{x) = I (0 1 d+ °=^ k) Vct=i( ^ eik ' x ' (4 ' 26b) 

and postulate that it annihilates any physical Fock state \^f) of the tachyonic field, 

(*\i>l-J 1 (x) = 4^ 1 (x)\*)=0. (4.27) 

These relations automatically imply that the Gupta-Bleuler condition also is realized in terms of the expectation 
value 



(4.28) 



but the condition (4.27) is stronger. We recall that the Gupta-Bleuler condition on the photon field reads 
(^■yld^A^l^^) = 0, where is a Fock state of the photon field. As stressed in Schweber's book [63] on p. 246, 
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the condition (^ r 7 |c? (tt A' 1 |\f r 7 ) = is not sufficient for the suppression of the longitudinal and scalar photons, but one 

must postulate that d^A^ |\& 7 ) = 0, where Ap is the positive-frequency component of the photon field operator. 
Because the tachyonic fermion, unlike the photon, is not equal to its own antiparticle, we need two conditions, given 
in Eq. (4.27). 

A crucial question now concerns the possibility of reversing the helicity-dependence, i.e., the question of whether 
or not a different choice for the helicity-dcpcndcnt factors in Eq. (4.3) exists that would imply negative norm for left- 
handed particles and right-handed ant ip articles. A related question is whether other tachyonic Dirac Hamiltonians 
exist for which the choice f(a, k) = g(a, k) = +<j instead of (—a) would fulfill our general postulate, namely, the sum 
rule (4.5). For reasons outlined in the following, we can ascertain that this is not the case; tachyonic spin-1/2 particles 
should always be left-handed. 

The arguments supporting this conclusion are as follows. First, if we assume that the tachyonic field fulfills a sum 
rule of type II which for massless fields is given in Eq. (2.6), then it is impossible to replace (— a) by (+cr) in the sum 
rule because of the necessity to preserve a smooth massless limit. The considerations in the text following Eq. (2.6) 
imply that if one were to replace (—a) by (+<r) in the massless case, then one would violate the sum rule (2.6). The 
second argument is obtained by explicit calculation. We have checked that if one replaces m — > — m in the tachyonic 
and imaginary-mass Dirac equations (1.7) and (1-9), then the sum rules fulfilled by the corresponding fundamental 
spinors still contain the characteristic factor (— a). For the imaginary-mass Dirac equation, this result is obtained in 
Ref. [25]. Intuitively, we can understand this result as follows: The mass m in the denominators of the right-hand 
sides of Eqs. (3.7) and (3.12) is obtained as the modulus \J m 2 — \m\ and does not change if we replace m — ¥ —m in 
the superluminal Dirac equation. The mass in the numerator of the right-hand sides of Eqs. (3.7) and (3.12) changes 
sign, but this is consistent with the obvious change in the functional form of the positive-energy and negative-energy 
projectors as we change the sign of the mass term. Again, this consideration supports the conclusion that we cannot 
invert the helicity-dcpcndcncc by choosing a different Hamiltonian; the factor (— a) persists. 

The third argument comes from the tachyonic Gordon identities discussed in Sec. IV C. We use the Gordon identity 
Eq. (4.23) and the normalization (3.6) to calculate the bispinor trace (with a 7 multiplied from the right) of the 
left-hand side of Eq. (3.7a), 



tr ]r(-a)^(fc)®W a (fc)7 5 7 - ^(- ( T)W ff (fc) 7 5 7 AI= °W (T (fc) 



= "E(- CT ) ( -—) U a {k)U a {t) = 2 - . (4.29a) 



The bispinor trace of the right-hand side of Eq. (3.7a) is 



^^) = <,(/AU 4 £ 2 * „ 29b) 



2m ) \ 2m J 2m m 

which shows the consistency of the bispinor sum (3.7a) with the Gordon decomposition (4.23). If we were to replace 
(— cr) by (+cr), the two sides of the relation (4.29a) would differ by a minus sign. The bispinor trace of the left-hand 
side of Eq. (3.7b) is 



tr (l> CT ) ® "Mfc) 7 5 7°) = X> CT ) Mk) ■ 



7 V= u V CT (fc) 

= E(- ff ) (-) V a {%) = 2 - . (4.30a) 
* — ' \ m J ■> v ' m 

— — (7 

We have used the tachyonic Gordon decomposition for negative-energy states as given in Eq. (4.23), which differs 
from the positive-energy Gordon decomposition by a minus sign, but an additional minus sign is obtained from the 
Lorentz-invariant normalization of the negative-energy fundamental bispinors. From the right-hand side of Eq. (3.7b), 
we have 

which again is fully consistent, but only because we have a (— a) in the sum rule (3.7), which combines with the 
(— a) from the Lorentz-invariant normalization of the V a {k)V cr {k), to give 2E/m as a final result in Eqs. (4.30a) 
and (4.30b). 
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V. PHYSICAL INTERPRETATION: FROM TACHYONIC NEUTRINOS TO COSMOLOGY 

A. Arguments for and against Tachyonic Neutrinos 

In the absence of conclusive experimental evidence, the hypothesis of tachyonic neutrinos has been controver- 
sially discussed in the literature. Three main arguments [64] have been brought forward against tachyonic neutrinos. 
(1.) They would require us to give up the notion of a Lorentz- invariant vacuum state, and even the vacuum would 

become unstable in the presence of tachyonic fields. (2.) Given the tachyonic dispersion relation E = \/k 2 — m 2 , the 
role of states with \k\ < m needs to be clarified. (3.) The physical (probability!?) interpretation of the conserved 
Noether current of the free tachyonic Dirac equation has been called into question [64] , and it has been argued that 
no consistent interpretation can be given because certain zero components of the conserved current were conjectured 
to vanish for all tachyonic momentum eigenstates [64] . 

A possible answer for question (1.) has been proposed in Ref. [9]. Summarizing the argument, it has been concluded 
in Ref. [9] that one can solve the problem in two ways, (i) One can Lorentz transform the vacuum state, and Lorentz 
transform all fundamental creation and annihilation operators of the fermion field (some of these will change from 
annihilators to creators upon transformation, due to the space-like nature of the tachyons) . (ii) One keeps a Lorentz- 
invariant vacuum state, and only transform the space-time arguments k^ and of the field operators, keeping all 
creation operators as creators and annihilation operators as annihilators. The amplitudes, cross sections, etc. obtained 
using approach (ii) then depend on scalar product of four-vectors which are equal to the result obtained by first 
calculating the process in the original Lorentz frame, and then, performing the Lorentz transformation into the 
moving frame. 

The conjecture regarding an expansion about a "false" vacuum in the presence of tachyons can be traced to the 
fact that most of the tachyonic theories discussed so far in the literature are scalar [28-33]. Indeed, scalar tachyons 
have a problem with instability, because of the structure of the mass term in relation to the field Hamiltonians, 
which changes sign m 2 — > —m 2 in a tachyonic theory, suggesting that the field energy can be lowered by creating 
tachyons. The problem does not occur in tachyonic spin-1/2 theories because a linear, not quadratic, mass term 
enters the field Lagrangian and Hamiltonian. Provided one reinterprets the spin-1/2 antiparticle solutions in the 
usual way (negative energy for propagation into the past becomes positive energy for propagation into the future), it 
then becomes immediately clear that the vacuum energy cannot be lowered upon spin-1/2 tachyon anti-tachyon pair 
production. 

An solution to problem (2.) has also been proposed in Ref. [9]. Namely, the energies with E = ±^/P — m 2 — i e 
(for \k\ < m) find a natural interpretation in terms of complex resonance and antiresonance energies, which describe 
unstable states which decay in time. Particle resonances are damped for propagation into the future, antiparticle an- 
tiresonances are damped for propagation into the past, as they should be [9]. The occurrence of momentum eigenstates 
with real energies, and resonances with complex resonance energies, is a well-known phenomenon all across physics 
(e.g., in atomic physics, the ground-state energy of the helium atom is strictly real, whereas auto-ionizing resonances 
in the three-particle system have a manifestly complex resonance energy). 

Another "myth" which should be refuted concerns a conceivable "runaway reaction" where a moving tachyon 
releases an arbitrarily large amount of energy, as it loses energy and accelerates, given its classical energy-velocity 
relation E — m/y/v 2 — 1. According to this relation, a tachyon indeed accelerates as its energy is lowered and becomes 
commensurate with the invariant mass square. (At high energy, a tachyon approaches the light cone, though.) An 
infinitely fast tachyon takes the role of a tardyon at rest [28] ; the explicit eigenstates have been indicated in Ref. [9] . 
However, energy conservation holds, and one cannot mix the notion of energy increase by acceleration, which only 
holds for tardyonic particles, in order to "convert" a tachyon losing energy by acceleration into a tardyon that gains 
energy in the same process. The "runaway reaction" is impossible; only a finite amount of energy is released as the 
tachyon accelerates and the energy goes from E = m/y/v 2 — 1 to zero. Energy conservation holds for tachyons, even 
if they accelerate, somewhat counterintuitively, when losing energy. 

An answer to question (3.) has not yet been provided so far in the literature to the best of our knowledge. Here, 
we aim to provide a possible physical interpretation for the conserved current and scalar product, and also show 
where certain arguments presented originally by Hughes and Stephenson in their research article [64] entitled "against 
tachyonic neutrinos" become inconsistent. The Lagrangian density of the tachyonic Dirac particle reads, in first 
quantization [6], 

C{x) = i>{x) 7 5 (i 7 M d M - 7 5 m) $(x) . (5.1) 
This is equivalent, up to partial integration, to the symmetric form [10] 



(5.2) 
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where we suppress the space-time argument x = (t, r). The non-symmetric form (5.1) clearly exhibits the presence of 
the "chiral adjoint" tp(x) 7 5 in the Lagrangian. The conserved current is 

J"(x) = #r) 7 5 -fiftx), d fl J»(x)=0. (5.3) 

The zero component 

J°(x) =^(x) 1 5 1 °^(x) (5.4) 
assumes the following values in plane- wave eigenstates, according to Eqs. (4.11a) and (4.11b), 

W CT (fc) 7 5 7°^W = -<?-, (5.5a) 

m 

V CT (fc)7 5 7°K(fc) = -a-. (5.5b) 

m 

We conclude that J° is positive for left-handed particle states {a = —1) and right-handed antiparticle states (likewise, 
a = — 1; the hclicity is equal to — a for antiparticles). For E — 0, the "axial norm" of the momentum eigenstates given 
in Eq. (5.5) vanishes, but this is a very special case. Namely, for tachyons, E — implies |fc| — m — mv/\/v 2 — 1 for 
the momentum and thus corresponds to an "infinitely fast" tachyon (v —¥ 00), which remains "infinitely fast" upon 
Lorentz transformation. It would thus be wrong to conclude, as done in Ref. [64] in the text following Eq. (54) of 
Ref. [64], that all tachyonic states have zero axial norm, because this would correspond to a forbidden generalization 
of a result which holds asymptotically, for v — > 00, to all finite values of the tachyonic velocity v. Our Eq. (5.5) gives 
the explicit result for any finite value of the energy E. 

The spatial integral of the zero component of the conserved current, 

d 3 rj°{x)= J d 3 r^(x)j 5 ^ip(x) = - J tP+(x)-/ 5 $(x) , (5.6) 

is precisely equal (up to a sign) to the scalar product (ipi(t),ip2(t)) = / d 3 r ipf(t, r) 7 5 ifaft, r) which is conserved 
under the time-evolution by the 7 5 Hermitian (pseudo-Hermitian) Hamiltonian H§ . This scalar product is not positive 
definite, as already noticed in the work of Pauli [13], and precisely corresponds to the scalar product introduced by 
Pauli in Eq. (3) of Ref. [13], where in the notation of Ref. [13] we have r\ = 7 s . Similar observations have been made in 
Eqs. (34) and (52) of Ref. [64]. According to Refs. [65, 66], one could otherwise define a so-called C operator, which is 
not equal to the charge conjugation operator, and "remedies" the problem of negative norm attained by some states 
under the 7 5 norm, leading to a redefined, positive-definite scalar product. However, the negative norm finds a rather 
natural interpretation, and a redefinition of the scalar product (tpi(t), ip2{t)) = J d 3 r tj)f(t, r) 7 5 fait, r) therefore is 
not required. 

For the ordinary Dirac equation, the conserved current is J fJ '(x) = ij) e {x) 7 M tp e (x). Its timelike component is 
J (x) — ip£(x)ip e (x), where the subscript e reminds us of the electron. The latter can be interpreted as a positive- 
definite probability density which is conserved under the time evolution generated by the ordinary Hermitian Dirac 
Hamiltonian H^K 

The tachyonic Dirac current J 11 is obtained from the ordinary Dirac by the replacement ip e ( x ) ~^ ipix)^ 5 . 
We have seen that the scalar product for the tachyonic Dirac Hamiltonian is equal to an integral over the timelike 
component of the conserved Noether current of the Dirac equation and is not positive-definite. In order to put this 
observation into perspective, it is instructive to recall that for the Klein-Gordon equation {d^d^ + m 2 )<p(x) — 0, the 
zero component of the conserved current j^{x) = (<p*(x)d fl (/)(x) — <p(x)d tl (j>*(x)) is not positive-definite, either. 
Therefore, the zero component of the Klein-Gordon current cannot be interpreted as a probability density but must 
be interpreted as a charge density, which is positive for particles and negative for antiparticles. 

This interpretation is not available for the zero component of the Noether current of the tachyonic Dirac equation, 
because the equation is not charge conjugation invariant and is primarily proposed to describe neutrinos [10]. However, 
we can come closer to a physical interpretation of the timelike component of the Noether current of the tachyonic 
equation if we compare the interaction Lagrangian £qed of quantum electrodynamics to the weak interaction Cw of 
a neutrino and a "heavy photon", i.e., a Z° boson, 

£qed= -e^VlM/i 6 ^(V^^W^ (5.7) 

Z Sin U\y COS U\Y \ z / 

where 6\y is the Weinberg angle. The axial vector part CA, of Cw is 

£ - = -2^M ?7V ^ Z - (5 - 8) 
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where the ordering of the 7 matrices is important. For the interaction with the time- like component of the vector 
potential, we have the expressions 

£°qed= -e?e7°i^o £&° = - 9 - L J 7 5 7V4 (5.9) 

For quantum electrodynamics (QED), we interpret ip e ~f°Tp e = ip+ tp e as the probability density, and this suggests an 
interpretation of the expression ip V' as an "axial probability density" or "axial interaction density" of the neutrino 
field with the timelike component of the Z° boson. According to Eqs. (5.5a) and (5.5b), the "axial interaction density" 
is positive for the physically allowed states (left-handed particle and right-handed antiparticle states), and negative 
for the physically forbidden states (right-handed particle and left-handed antiparticle states). This consideration is 
independent of the suppression mechanism for the states of "wrong" helicity which, in second quantization, due to 
negative norm [see the discussion following Eq. (4.25)]. 

In general, the physical interpretation of tachyonic theories has been discussed by Feinberg [32, 33] and Sudarshan 
et al. [28-31]. The rcintcrpretation principle is a cornerstone of the theory. The only physically sensible quantities in a 
quantum theory are transition amplitudes. If the time-ordering of superluminal events changes upon Lorentz transfor- 
mation, then one reinterprets the amplitude as connecting two space-time events whose coordinates are transformed 
according to the Lorentz transformation, so that the only physically sensible quantity (transition amplitude) simply 
connects two events one of which happens before the other [31]. This may seem counter-intuitive at first, but it is 
perhaps a little less counter-intuitive if we take into account that the accepted formulation of quantum field theory 
is based on the Feynman propagator and on the reinterpretation principle for the "advanced" part of the Feynman 
Green function, which propagates anti-particle solutions into the past. In order to avoid problems with regard to 
causality, the canonical quantum field theory of subluminal particles has to be supplemented by a reinterpretation 
principle, just like the tachyonic theory. The physical interpretation of the Noether current is not affected by this 
consideration. 

Finally, let us point out that even without invoking reinterpretation, superluminal propagation can be compatible 
with causality if we postulate that the tachyonic mass m is so small that the superluminality is within the limits set 
forth by the uncertainty relation. With an energy E — m/ f3 2 — 1 w m/V26 with v = 1 + S and 5 <C 1, we have 
AE At w m/V26 At < 1 (in units with h = 1). Thus, for an infinitesimal tachyonic mass parameter which does not 
exceed m < y/25/At, causality is preserved within the limits set by the uncertainty principle. Quantum limitations, 
the role of unstable modes and quantum tunneling in superluminal propagation have been discussed in Refs. [67-70]. 



B. Gupta— Bleuler Condition and Seesaw Mechanism 



The commonly accepted mechanism for the suppression of right-handed neutrino and left-handed antineutrino 
states is the seesaw mechanism [71]. After integrating the heavy degrees of freedom (the sterile neutrino), it contains 
a nonrenormalizable dimension-five operator and has a hierarchy problem: Namely, the neutrino masses are inversely 
proportional to the grand unification (GUT) scale A and sensitively depend on fine-tuning of A. Let us consider 
Eq. (18) of Ref. [71] and hypothetically consider a formulation that would result if the physically observable neutrino 
were right-handed. Then, we could reformulate Eq. (18) of Ref. [71] as 



Cf = - 



Rt'rH 



H T (RtrY 



h.c, 



(5.10) 



where 
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With the expectation value of the Higgs field, 



J_ ( v 

instead of the left-handed Majorana neutrino masses, the right-handed ones would be small, 

£M = - \ J2 9 r>RM? r , {v rR f + h.c. , M£, = ^ Y r , r . 



(5.11) 



(5.12) 



(5.13) 
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The seesaw mechanism is not unique in suppressing a definite helicity of the neutrino. It is unique provided one 
formulates it in terms of the left-handed fermion fields, but it could be formulated with inverted helicities if the 
helicity of the observed neutrinos at low energy were different. In the latter case, the right-handed neutrino mass 
would be small, instead of the left-handed one. On the other hand, the seesaw mechanism has the distinct advantage 
that it is not necessary to assume a superluminal character of the neutrino. 

The mechanism discussed here in the text following Eq. (4.25) is definite in making a prediction regarding the 
suppression of right-handed neutrino states, as explained by three independent arguments in Sec. IV D. However, we 
have to assume a superluminal neutrino. While an interacting superluminal field theory is problematic, it is perhaps 
not as problematic as previously thought (see Ref. [72] and Sec. 4 of Ref. [9]). 

Final clarification can only come from experiment. The seesaw mechanism is compatible with a Majorana neutrino. 
The tachyonic Dirac equation implies that the neutrino cannot be equal to its antiparticle, because it does not allow 
charge-conjugation invariant solutions. It is only CV, but not C invariant. Experimental evidence for neutrinoless double 
beta decay is disputed [73] , and direct measurements of the neutrino mass square currently exclude neither positive 
nor negative values [44-50]. The generally accepted seesaw mechanism implies that neutrino masses are generated 
by a nonrenormalizable interaction with a concomitant hierarchy problem, and the mechanism in itself could be 
reformulated with opposite helicities. It is compatible with a Majorana neutrino. By contrast, a tachyonic neutrino is 
"automatically" left-handed and not equal to its own antiparticle. It is described by a 7 s Hermitian Hamiltonian and 
is plagued with the conceptual difficulties associated with (ever so slightly) superluminal propagation. This means 
that it is experimentally possible to test the models. 



C. Tachyonic Neutrinos as a Candidate for Dark Energy 

The formulation of a gravitational interaction of a spin-1/2 particle is nontrivial in the quantized formalism. Brill 
and Wheeler [74] performed the pioneering steps in this direction. In order to formulate the gravitational coupling of a 
Dirac particle, one has to formulate generalized Dirac matrices 7^, which fulfill anticommutation relations compatible 
with the local metric g^" (x) of curved space-time. Based on the Christoffel symbols T p ^ = T p M „(x), one formulates the 
Christoffel affine connection matrices T M in spinor space and calculates the covariant derivative as follows [74-78], 

{fW.fWl = 2 <T(z) , rv = \f (§^f + ^ - |^) , (5.14) 

v^ = ^-r^%+%r,-r,7 M = o, ^ = -\? f^-^r^) . (5.15) 



The above formula for is valid in the case of a diagonal metric ~g^ v such as the Schwarzschild metric con- 
sidered in Ref. [76]. For a general space-time metric, and with "West-Coast" conventions for the local vierbein 
g»" = diag(l, — 1, —1, —1) ("East-Coast" conventions were used in Ref. [74]), the result reads as 

r *=-^ (^ a ^- r%fc )^' (5 - i6) 

where = i [7^,7"] is the spin matrix. The a a p and bj^ coefficients transform the Dirac matrices to the local 
vierbein, 

7„ = V"7a, 7p = «%7 q , 7 Q =«%7 P , l a = b p a i p - (5.17) 
With the covariant derivative V M = <9 M — r„, the gravitationally coupled Dirac equation reads as [74-76] 

(i7"V M -m)ip{x) = 0. (5.18) 
In the case of a tachyonic Dirac particle, it has to be reformulated as follows, 



(i7 M V M -7 5 (x)m) V(a:) =0. 



(5.19) 



The space-time coordinate-dependent matrix 7 (x) can be defined as 



_ 5(x) = ^^ r _ VT5= l_ - a -P-,-S = if >-l-2-3 (5 2Q) 
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where g = det is the determinant of the metric, and e Q ^ 7 5 = £ Q ^ 7 <5 is the local e tensor, while £ Q /3 7 5 is the 
totaly antisymmetric Levi-Civita tensor. The last identity in Eq. (5.20) is valid for a diagonal metric ■ 

Let us first discuss Eq. (5.18) very briefly. A projection onto the upper and lower radial components /(r) and g(r) in 
a gravitational field can be found in Eqs. (19) and (20) of Rcf. [76]. For vanishing electrostatic potential V — > 0, Eq. (20) 
of Rcf. [76] is invariant under the replacement f(r) <H> g(r), and E <-> — E. So, if E is an eigenvalue of the gravitationally 
coupled Dirac equation, so is — E. Invoking reinterpretation and replacing — E — > E for antiparticles, we find that 
the spectrum of the gravitationally coupled Dirac Hamiltonian is the same for particles and antiparticles. Therefore, 
the formalism makes the unique prediction that tardyonic antiparticles, like tardyonic particles, are attracted by a 
gravitational field. In passing, we note that the often cited motivation for the investigation of trapped antihydrogen 
and its interaction with the gravitational field therefore is faced with a unique theoretical prediction: Antiparticles 
are attracted by gravitation as much as particles are. 

It has been confirmed within the last two decades [79-82] that the Universe expands more rapidly on large distance 
scales than compatible with the matter density in the Universe. Coupling to a scalar field ("quintessence") is usually 
invoked in order to explain the expansion of the expansion rate of the Universe [83] . As the scalar quintessence field 
"rolls down its potential" , it accelerates the expansion rate of the Universe. Because of a self-attracting property (the 
quintessence field energy can be lowered by increasing the local density of the quintessence field), quintessence has 
positive energy density but negative pressure. This property is necessary in order to constitute a candidate for dark 
energy [83]. In the standard model of cosmology, dark energy accounts for about 73% of the total mass-energy of the 
Universe. The quintessence field thereby acts like a time-dependent cosmological constant. 

We can thus conclude that most of the energy in the Universe actually is not gravitationally attractive, i.e., that 
gravity can repel. In the following, we shall present qualitative arguments which suggest that tachyonic neutrinos 
may play a role in the expansion of the Universe, and may contribute to "dark energy" . Conceivable connections of 
tachyonic physics and dark energy have been explored in the literature, employing either scalar fields [84] or Lorentz- 
violating mechanisms [85]. In order to provide an alternative explanation for dark energy, it is necessary to invoke a 
mechanism that leads to a repulsive gravitational force on intergalactic distance scales in the Universe. It is not fully 
surprising that tachyonic neutrinos may provide for such an alternative mechanism. Namely, both on the classical 
level [86], as well as on the level of quantum theory [Eq. (5.19)], tachyonic particles are repulsed by gravitational 
fields. 

On the classical level, this is seen as follows. We start from the familiar equation of motion of a particle with mass 
to in curved space-time [see Eq. (86') of Ref. [86]], which is a geodesic, 

d 2 ^ „ dx p dx° 

^7 + r ^^^7 = °- (5 - 21) 

With a suitably redefined proper time ds', the zero geodesic for a tachyon reads as [see Eq. (86') of Ref. [86]] 
The force exerted on a tardyon reads, according to Eq. (79a) of Ref. [86], 



d 2 r^ dr p dr a 

F * = m Vir = - mT %°^r^-- ( 5 - 23 ) 

ds z ds as 



However, there is a sign change for a tachyon [see Eq. (79b) of Ref. [86]], 

d 2 x'P _ mr , M dx'P c 
ds' 2 P<T ds ds 



d 2 r' M dx'P dr" 7 

F " = -™ %V = ™ r V -r- . (5-24) 



corresponding to the change in the energy- velocity relation E = m/yl — v 2 — > E = rn/\/v 2 — 1 for tardyon ver- 
sus tachyon, which leads to gravitational repulsion for tachyons. Another intuitive way to understand gravitational 

tachyonic repulsion is to observe that the quantity E = Vk 2 + m 2 gets smaller when the position-dependent mass 

to — m(x) gets smaller, whereas the tachyonic energy E = \/k 2 — m 2 gets larger when the position-dependent mass 
to = m(x) gets smaller. 

As neutrinos are ubiquitous within the cosmos, one now needs to evaluate their conceivable contribution to the 
repulsive force on intergalactic distance scales. A detailed discussion is beyond the scope of the current article, but we 
can formulate some initial considerations and order-of-magnitude estimates. One can identify, a priori, two sources of 
neutrinos which may become important, (i) thermalized neutrinos which decouple in the early Universe [87-90], and 
(ii) non-thermalized, high-energy neutrinos which are continuously generated through various cosmic processes, such 
as nuclear fusion in stars like the sun, supernovae, and high-energy cosmic background. 
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It is generally assumed that today's cosmic background (CMB) radiation is accompanied by a neutrino background, 
which at some point (low energy) decoupled from the other particles. The average energy of the background neutrinos 
is generally assumed to be of the same order-of-magnitude as compared to the cosmic microwave background (CMB), 
and the currently accepted value for the temperature of the neutrino background is (4/11) 1 / 3 = 0.714 times that 
of the electromagnetic background radiation, i.e., about 1.9 K. The Universe, according to this hypothesis, would be 
filled with a sea of nonrelativistic background neutrinos. This picture changes when we assume that the neutrino 
is tachyonic. In the earliest stages of the Universe, tachyonic neutrinos are of course highly energetic, with E = 
\/k 2 — to 2 w \k\ ^> 0. As they lose energy, they become faster, because E = m/\/v 2 — 1 — > (in the classical 

theory). The energy E = VW — to 2 vanishes as |fc| — > to, within quantum theory. For \k\ < to, according to [9], the 
spectrum of the tachyonic Hamiltonian contains unstable resonances and anti-resonances, with a purely imaginary 
resonance energy. Because the real part of the resonance energy E = — iy m 2 — k 2 vanishes, the unstable states form a 
"background" of fluctuating quantum states in the Universe. Particle resonances (E — —\\J m 2 — fc 2 ) and antiparticle 

antiresonances {E = +iv to 2 — fc 2 ) are damped for propagation into the future and past, respectively (see Ref. [9]). 
The latter case amounts to propagation into the future if one invokes reinterpretation for the antiparticle solutions [9] . 

A very course-grained estimate regarding the role of the fluctuating resonance and antiresonance states can be given 
as follows. We start from the grand canonical emsemble 



n = -2^ V J q dfc k 2 in ^1 + exp ^^-^ ) j , (5.25) 

where ks is the Boltzmann constant, V is the normalization volume, p is the chemical potential, and E(k) is the 
energy of a tachyon as a function of k = \k\. The multiplicity factor 2 takes into account the particle and antiparticle 
resonances and is supplemented here in comparison to the derivation presented in Ref. [91]. Furthermore, we have 
assumed an isotropic energy dependence J d 3 fc — > 4ir / dfcfc 2 . The pressure is found as 

,-n-^f«^(i + «p(^)). 

For T — > 0, the contribution of resonances and antiresonances with |fc| < to to the pressure and energy density is 
easily evaluated according to Eqs. (27) and (28) of Ref. [91]. For the pressure, we have 

i f^odlmffi) i f m „ k 4 Ito 4 

Po = ^r dkk 3 — ^ = — - / dfc - , 5.27 

3tt 2 J dk 3ir 2 J Q ^ m 2 - k 2 16tt V ' 

whereas the energy density is 

Po = -^ dkk 2 lmE(k) = / dfc k 2 Vm 2 - k 2 = = -p . (5.28) 

7T 2 J Q TT 2 J 16 7T 

The equation of state fulfilled by the zero-temperature limit of the tachyonic Dirac sea is wq — po/ po = — 1. This is the 
required equation of state for a "vacuum" energy density that describes a nonvanishing cosmological constant [92] . If 
Po + Po = 0, then there is no net energy gain upon pulling on a "piston" which contains the Universe (an illustrative 
discussion can be found in Ref. [93]). 

Both results in Eqs. (5.27) and (5.28) are imaginary. This raises the pertinent question of how to incorporate the 
imaginary pressure and energy density of the neutrino field into the evolution equations of the Universe. In giving 
rise to unstable resonance states, the tachyonic Dirac field, if it exists, would be different from any other known 
fundamental quantum fields [59] . An exhaustive mathematical description would require an extension of scattering 
theory to complex energy and momentum exchanges, applied to low-energy collisions with neutrinos, and describe the 
continuous decay and repopulation of the unstable tachyonic states in the Universe. A detailed discussion is beyond 
the scope of the current article. However, in a first approximation, we can argue that an established technique for the 
mathematical treatment of complex resonances entails complex scaling (for reviews, see Refs. [94, 95]), and we thus 
explore a complex scaling ansatz to the solution of the evolution equations in the following. 

The imaginary results in Eqs. (5.27) and (5.28) describe the quantum fluctuations due to the unstable states. We 
thus consider Eqs. (8)— (11) of Ref. [92], which relate the accelerated expansion rate a/\a\ and the Hubble constant 
H = a/\a\ (the dot denotes differentiation with respect to time) to the energy densities of various cosmologically 
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relevant quantities, 



8ttG , , 8ttG 



a 87rG 4-7rG . 47rG , /in 

rr = — ^— PA ^— (Pm + 3pM + Pk + apt) ~ —z— (2pA - PM) ■ (5.29b) 

\a\ 3 6 3 

Here, p is the matter density in the Universe, pk — — 3fc/(87rGa 2 ) is the energy density associated with the curvature 
of the Universe [k = +1,0,-1), while pa is the energy density corresponding to the cosmological constant. The 
mass density is pm and pm ~ 0. We aim to solve Eq. (5.29) using a complex variable ansatz a = \a\ exp(i^) and 
t = t exp(vp). Using Eq. (5.29), with a value H — 71km s _1 Mpc -1 for the Hubble constant, one easily reproduces 
the accepted value for the critical mass density of p cr i t ~ 9.34 x 10~ 33 kg cm -3 (in SI units) which corresponds to 
a critical energy density of p cr i t « 8.39 x 10~ 9 J m~ 3 (we have set c = 1). In natural units, the result converts to 
p crit « 9.99 x 10- 45 GeV 4 . 

We courageously set k = in Eq. (5.29) (flat Universe) and match the modulus of the (complex) accelerated 
expansion a/\a\ against the currently accepted values of Am = p/pcrit ~ 0.27 (see Ref. [96]) and — Pa/Pci-h ~ 0.73. 
According to Eq. (5.29b), the matching is performed by courageously assuming that the bulk of the cosmological 
constant is given by the complex energy density (5.28), and equating the complex modulus as follows, 



PA-Po^-- 7r , \2p' A - PM \ = \J4\p'a\ 2 +Pm = 2/j A -pM«l-19pcrif (5.30) 

ID 7T v 

p' A = — i 0.579 Pcrit, and so our estimate for the neutrino mass reads as 



(5.79 x 10- 4b GeV 4 ) - \p \ = — m~ 0.0232 eV. (5.31) 

1D7T 1 1 



If this treatment is valid, then the mass m necessarily has to be the mass of the "heaviest" neutrino mass eigenstate, 
i.e., the one with the largest modulus of the tachyonic mass m. Heavier eigenstates would lead to a larger cosmological 
constant, because po is proportional to the fourth power of the neutrino mass. In terms of a conversion to flavor eigen- 
states, the heaviest mass eigenstates might well be the one closest to the electron neutrino and electron antineutrino 
lepton flavor eigenstates [97]. 

We have interpreted the quantities po and po given in Eqs. (5.27) and (5.28) as contributions to a cosmological 
constant which describes the evolution of the Universe in complex space and time directions a = |a| exp(i6>) and 
t = t exp(iys). We should supplement the solutions for the complex rotation angles for space and time, which read as 
9 = 38.1° and cp = 70.6°, within our complex scaling approach. The evolution of the energy density with the aging of 
the Universe is given by [see Eq. (12) of Ref. [92]] 

po = pi(0)a- 3{1+wo) = const., w = -l. (5.32) 

We therefore obtain a time- independent energy density po, which is consistent with the fundamental character of the 
fluctuating unstable resonances of the neutrino field. We reemphasize that the mass m ~ 0.0232 eV is tachyonic and 
enters the gravitationally coupled tachyonic Dirac equation (5.19). The given mass value is not currently excluded by 
any terrestrial experiment, and it is consistent with the observed time spread of the arrival times of neutrinos from the 
supernova SN1987A (see Ref. [98]). The hypothesis that very light elementary spin-1/2 particles could be tachyonic 
has been pursued elsewhere [99]. 



VI. CONCLUSIONS 



In the current investigation, we present the fundamental solutions of generalized Dirac equations in the helicity 
basis, in a systematic and unified manner. Of particular importance are the Dirac equation with two tardyonic mass 
terms mi -M7 5 1TI2 and two tachyonic mass terms imi +7 5 rri2. Let us summarize the main results. We have discussed 
the ordinary (tardyonic) Dirac equation in Sec. II B, a tardyonic Dirac equation (with two mass terms) in Sec. II C, 
and two tachyonic Dirac equations in Sees. Ill A and IIIB. We give the fundamental eigenspinors that enter the 
plane-wave solutions of all of these equations [see Eqs. (2.10), (2.11), (2.17), (2.18), (3.2), (3.3), (3.9) and (3.10)]. 
The eigenspinors are obtained using projector techniques as outlined in Chap. 2 of Ref. [59]. For the "normal" Dirac 
equation [see Eqs. (2.10) and (2.11)], our results are consistent with Ref. [3] and Chap. 23 of Ref. [60]. For the 
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generalized Dirac equations, the solutions have not yet appeared in the literature in the compact form given in the 
current article, to the best of our knowledge. For the ordinary (tardyonic) Dirac equation, for the tachyonic Dirac 
equation, and for the imaginary-mass Dirac equation, the prefactors are brought into compact analytic form [see 
Eqs. (2.10), (2.11), (3.2), (3.3), (3.9) and (3.10)]. 

Finally, in Sees. IV A and IV B, we hnd that the tardyonic and tachyonic theories can be unified on the basis of the 
structurally simple anticommutator relations given in Eqs. (4.14) and (4.15), which are independent of the magnitude 
of the mass terms. As outlined in Sec. IV A, the coefficient functions / = f(cr,k) and g = g(a,k) in the postulated 
form of the anticommutators (4.3) enter the tensor sums over the fundamental spinors in Eq. (4.5). For the egregiously 
simple choices indicated in Eqs. (4.14) and (4.15), which are consistent with a smooth massless limit, the tensor sums 
over the fundamental cigenspinors yield the positive-energy and negative-energy projectors, for both tardyonic as well 
as tachyonic eigenspinors. Consistency with the massless limit requires the presence of the factor (— a) in the tensor 
sums over the eigenspinors for the tachyonic equations. This fact is verified, on the basis of tachyonic Gordon identities 
and related considerations, in Sees. IV C and IV D. The presence of the factor (—a) in the fundamental tachyonic field 
anticommutators in Eq. (4.15) implies the suppression of right-handed particle and left-handed antiparticle states, 
due to negative norm, as shown in Eq. (4.25). 

Finally, in Sec. V, we observe that since tachyons are repelled by gravity, it might be worthwhile to investigate 
their conceivable role in the mechanism(s) responsible for the accelerated expansion of the Universe ("dark energy"). 
Furthermore, the tachyonic resonances and anti-resonance energies might play a role in the sum over states that enters 
the thermodynamic potentials of a free tachyonic fermionic gas in the low-temperature limit. If we consider the states 
with imaginary energy to be unstable, fluctuating states, then it is intuitively obvious that they might contribute to a 
fluctuating energy density and pressure on large distance scales in the Universe. As described in a somewhat approach 
in Sec. VC, an order-of-magnitude calculation based on a complex scaling transformation of the time evolution of 
the Universe leads to a tachyonic neutrino mass which is not excluded at present by any terrestrial experiment. We 
believe that it might be worthwhile to explore the physical consequences of the "ugly duckling" (tachyonic neutrino) 
somewhat further. It allows us to retain, among other things, lepton number conservation as a symmetry of nature. 
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What happens if Dr. Spock overtakes a left-handed neutrino.. 




Spinor Sum for Positive-Energy States: 
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Spinor Sum for Negative-Energy States: 
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There is a paradox which results if one overtakes a left-handed neutrino, looks back and sees 
the same particle right-handed, because right-handed neutrinos have never been observed in 
nature. One possible way to solve the paradox is sketched in the current article; the spinor 
sums are given in Eq. (3.7). Further explanations are in the text. 



